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Abstract 

It has been argued recently that parton showers based on colour dipoles conflict with 
collinear factorization and do not lead to the correct DGLAP equation. We show that this 
conclusion is based on an inappropriate assumption, namely the choice of the gluon energy 
as evolution variable. We further show numerically that Monte Carlo programs based on 
dipole showers with "infrared sensible" evolution variables reproduce the DGLAP equation 
both in asymptotic form as well as in comparison to the leading behaviour of second-order 
QCD matrix elements. 



1 Introduction 



Due to the large centre-of-mass (CM) energies at modem hadron colliders, and especially in 
connection with the production of high-f states it is becoming increasingly essential that we have 
well-controlled calculations that cover all regions of phase space in a reliable way (see, e.g., [1]), 
and that these are implemented consistently in the phenomenological tools we use for collider 
studies — we shall here focus on Monte Carlo Event Generators. It is likewise important that 
we have an accurate theoretical understanding of the properties we rely on these calculations to 
have. 

Non-trivial collider observables almost invariably involve an interplay between widely sep- 
arated energy scales. This represents an important challenge to collider phenomenology since, 
in any gauge theory with massless gauge bosons (thus QED and QCD in particular), such scale 
hierarchies give rise to logarithmic enhancements in the matrix elements, order by order in per- 
turbation theory, which ultimately render a truncation of the perturbative series invalid at any 
fixed order. This becomes increasingly relevant as the collider CM energy (or other relevant hard 
scale) grows, leaving more room for scale hierarchies to develop beneath it. 

A state-of-the-art collider physics calculation includes both a good description of physics 
at short distances — usually represented by perturbative leading-order (LO) or next-to-leading 
order (NLO) matrix elements with renormalization-group (RGB) improved couplings — a good 
description of the transition from short to long distances, which takes proper account of any 
large scale hierarchies that may develop on the way — usually represented by leading-log (LL) 
parton showers incorporating as many next-to-leading-log (NLL) effects as feasible, as well as 
models of other possible perturbatively enhanced aspects such as multiple parton interactions 
(MPI) — and finally a good description of the physics at long distances — usually represented 
by non-perturbative models of beam remnants, hadronization, and hadron decays. 

We shall here focus on parton shower algorithms, which provide the connection between 
the perturbative fixed-order matrix elements and the non-perturbative hadronization models, 
and which thus constitute an essential ingredient of general-purpose event generators like Her- 
WIG [2,3], Pythia [4,5], or Sherpa [6]. 

Parton showers generate infinite-order approximations to matrix elements (both real and vir- 
tual), in such a way as to coincide exactly with the matrix elements in the singular limits. The 
number of singular coefficients that are reproduced exactly depends on the order of the par- 
ton shower; thus, an LL parton shower can be expected to generate the correct coefficients for 
the "leading" matrix-element singularities (to be defined further below). These terms dominate 
in the limit of infinitely large hierarchy between the scales of each successive emission, and 
hence LL showers are supposed to be exact in this so-called "strongly ordered" limit. An NLL 
shower [7-10] should also generate the correct coefficients for the next-to-leading singular terms 
(dominant in regions with one less large hierarchy than the strongly ordered limit), and so on. 
While the accuracy of the above mentioned general-purpose event generators can be debated, 
we note that none of them systematically include the tree-level n ^ n + 2 and 1-loop n ^ n + 1 
splitting functions, and hence they are all formally LL. However, also without exception, they do 
incorporate a number of non-trivial NLL effects systematically and usually perform significantly 
better than corresponding analytical LL calculations. Any proposed algorithm should therefore 
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be subjected to two basic tests: 1 ) whether it correctly reproduces QCD in the LL singular limit, 
and 2) how well it approximates QCD in NLL singular limits. 

At first order in perturbation theory the matrix elements become singular in phase space re- 
gions corresponding to the emission of coUinear or soft particles. The first showering algorithms 
started from the coUinear factorisation of the matrix elements and approximated colour inter- 
ference effects through angular ordering [11-13]. An alternative approach is the Lund-dipole 
or dipole- antenna shower model, first implemented in ARIADNE [14-18], in which first-order 
colour interference effects are instead taken into account by choosing colour-connected pairs of 
partons to be the radiating entities. 

Here, unfortunately, a digression on nomenclature is necessary. These parton pairs were sim- 
ply called dipoles in the original Lund papers [14-18]. They are also known as antennae [19-21], 
mostly in fixed-order contexts. In order to be clear about what we mean to both communities, 
we therefore use the term dipole-antennae [22] here for these objects. Later, Catani and Sey- 
mour invented a related, but different, concept which they unfortunately also called dipoles [23]. 
Roughly, a Catani-Seymour dipole corresponds to one half of a dipole-antenna, partitioned in 
a way that isolates the coUinear singularities; hence the Catani-Seymour dipole could also be 
called a partitioned-dipole, as suggested in [24]. The problem has been further compounded by 
a new generation of shower models based on Catani-Seymour dipoles being generally referred to 
as "dipole showers" [25-28], whereas the same name was already used by Ariadne to describe 
its final-state radiation (FSR) model based on dipole-antennae [18]. To make the confusion com- 
plete, while Ariadne's algorithm for final-state radiation is based on a strict dipole-antenna 
picture, its treatment of initial-state radiation (ISR) is somewhat different. We shall nonetheless 
attempt to retain some measure of clarity by referring consistently to the Catani-Seymour type 
as partitioned-dipoles and to the antenna-type as dipole-antennae. The two shower types do ex- 
hibit many similarities, so the distinction may mostly be important to experts, but as with all 
approximations, the devil is in the details. 

The last few years have witnessed significant progress in the improvement of parton shower 
algorithms. Based on a proposal by Nagy and Soper [25,26] new algorithms based on partitioned- 
dipoles have been developed [27, 28]. A new final-state algorithm relying on the dipole-antenna 
picture has also been constructed [22], with similar properties as the Ariadne final-state shower, 
and a complete dipole-antenna shower model incorporating both ISR and FSR has been devel- 
oped in [29] . These new algorithms have their origin in subtraction methods for fixed-order cal- 
culations [23, 30-34] and implement in a correct way simultaneously the soft and the coUinear 
limit in a manner similar to that of Ariadne's FSR model. At the same time the new partitioned- 
dipole algorithms (as well as any dipole-antenna algorithm) are able to satisfy simultaneously at 
each step momentum conservation and the on-shell conditions. This is possible, because they are 
based on 2 ^ 3 splittings, where the spectator can absorb the recoil. Within the traditional 1^2 
splitting algorithms it is impossible to satisfy simultaneously momentum conservation and the 
on-shell conditions for a splitting. What has to be done in 1 ^ 2 splitting algorithms is either to 
restore momentum conservation by an ad-hoc procedure in the end, as in Herwig, or to restore 
it more locally during the shower evolution, but still involving at least a third parton at each step, 
as in Pythia and Sherpa. 

In a recent paper Dokshitzer and Marchesini [35] study a soft dipole-type model with recoil 
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effects, which is obtained from considering multiple antennas in QCD [36]. The first version of 
ref . [35] on the archive claimed that dipole showers are in conflict with coUinear factorization and 
do not lead to the correct DGLAP equation. Ref. [35] motivated us to study this question in detail. 
In this paper we show that the above mentioned conclusion is a consequence of an inappropriate 
assumption, namely the choice of the gluon energy as evolution variable. We further show that 
with a proper evolution variable the DGLAP equation is reproduced, thus proving that dipole 
showers, whether of the partitioned or the antenna type, have the correct LL behaviour, so long 
as sensible evolution variables are chosen. Subleading differences between these shower models 
will still be present, at the NLL level, which can first be accessed at 2"'^ order in perturbative 
QCD. We therefore also include a set of explicit comparisons of different shower models to 2"'^ 
order QCD matrix elements. We note that in a recent paper Nagy and Soper [37] have given 
a strict formal derivation that dipole showers reduce to the DGLAP equation in the strongly 
ordered limit. We further note that the second version of ref. [35] on the archive only claims that 
models with the gluon energy as the evolution variable conflict with coUinear factorization. This 
is in line with the findings of ref. [37] and of this paper. 

This paper is organised as follows: In sect.|2]we discuss the factorization of tree-level matrix 
elements in soft and coUinear limits. In sect. |3] we review the way a parton shower is obtained 
from the factorization properties of the matrix elements. We discuss in detail the choice of 
the evolution variable and point out that an energy-ordered shower is not compatible with the 
coUinear limit. Showers ordered by the transverse momentum or the virtuality are unproblematic. 
In sect. |4] we consider the evolution of the non-singlet quark fragmentation function. In sect. [5] 
we compare the analytical result of the previous section with numerical results obtained from 
Monte-Carlo programs based on both parti tioned-dipole and dipole-antenna showers. Again we 
show that, so long as "infrared sensible" evolution variables are chosen, these showers correctly 
reproduce the coUinear limit. Finally, sect. |6] contains our conclusions. 



2 Basics 

To set the scene let us consider the matrix element squared for y* — q{pi)g{p2)q{P2) in four 
dimensions: 

l-^3(Pi,P2,P3)| = ^egNcCF[2 H H . (1) 

V ■^12*23 523 SnJ 

In this formula, e denotes the electro-magnetic coupling, g the strong coupling, Nc = 3 the 
number of colours and Cf = {N^ — l)/{2Nc). The invariants are Sij = {pi + PjY and stjk — 
{Pi + Pj+Pkf- 

In the limit where the momentum p2 of the gluon becomes soft, the formula factorises as 

lim |j?3(pi,p2,;?3)P = '&nasCF^iKpi,P2,P3) \^2{p\,P3)\^ : (2) 

where 

B\k{pi,p2,p?,) = 2^^, (3) 

*12*23 
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and tti- = I (4;r). The matrix element squared for y* — > q{p\)q{p'i) is given by 

\^i{p\,Pi)'^ = 4e^NcSi2- (4) 

In the limit where the momentum p2 of the gluon becomes coUinear with the momentum pi of 
the quark such that p\=zP and p2 = (1 — z)^ we have the factorization 

Y\m\^3{pi,P2,P3)\^ = ^nasCF Pq-^qg\^2{P,P3)\^ , (5) 

Pl||P2 

with the Altarelli-Parisi splitting function 

A similar factorization formula holds for the case where the gluon becomes coUinear with the 
antiquark. 

With the help of an antenna function [21, 30] we may combine the three singular limits {p2 
soft, P2\\pi and P2\\P3) into one formula: 

lim \ji3{pi,p2,P3)\^ = SKasCfA^{pi,p2,P3) 1^2(^1,^3)^- (V) 
P2 unresolved 

The antenna function is given by 

A3{phP2,P3) = 2 H + , (8) 

*123 V *12'^23 ^^23 ■^12/ 

which is exactly the object used in both Ariadne and default Vincia for qq qgq branchings. 
The momenta pi and p3 entering the matrix element J^2 are obtained from p\, P2 and p^ such 
that they approach the correct limit in all singular limits, one possibility is [30]: 

(l+p)5l23-2r523 , , (1 -p)5l23-2r5i2 

Pi = ^7 ^ Pl+rp2^ r P3, 

2(5123-^23) 2(^123-512) 

(l-p)5l23-2(l-r)523 . , (l+p)5i23-2(l-r)5i2 

P3 = ^7 ^ Pl + {l-np2^ —, ^ P3, (9) 

2(5123-523) 2(5123 -512) 

where 



p = Jl+Ml-r)^. (10) 

512 +523 V 5123513 

Note that eq. O gives the correct factorization in any singular limit: p2 soft, p2\ \pi and P2\\P3- 
Let us introduce the dimensionless quantities 

512 1 523 1 ■^IS 1 /I IX 

yi2 = = 1-^3, 3^23 = = l-Xi, yi3 = = 1-JC2, (H) 

5123 ^^123 ^^123 
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Figure 1 : The Dalitz plot for the phase space of the unresolved particle. The location of the soft 
singularity is shown in red, the location of the coUinear singularities in blue. 

where the Xi are the ordinary energy fractions, evaluated in the CM of the 3-parton system, 

Xi = . (12) 

Obviously we have yn+yis +yi3 = 1 and X1+X2+X3 = 2. The antenna function reads then 

^3(Pl'^2,P3) = 2 H + . (13) 

■^123 V 3'12y23 3^23 ^12/ 

The unresolved phase space in terms of these variables is 

d<^unresolved = Y^^^3'12^^}'23^^}'138 ( 1 - ^12 + J23 + Jb) ) 0(yi2)0(j23)0(jl3) ■ (14) 

The Dalitz plot for the unresolved phase space is shown in fig [TJ The allowed phase space con- 
sists of the triangle yi2 > 0, 3^23 > 0, yn + J23 < 1- The soft singularity corresponds to the point 
yn = yis =0, while the coUinear singularity pi\\p2 corresponds to the line yn = (and the 
coUinear singularity piWps corresponds to the line 3^23 = 0). 

We can write the antenna function A3 in terms of two dipoles: 

A'^^{pup2,P3) = Dn3+D32,u (15) 

with 

(^^^+-). P3.,^^( y -^). (16) 

*i23 \yi2{yi2+y32) ynj sm Kyisiyn+yis) yn) 

In their recent paper Dokshitzer and Marchesini [35] study a model consisting of only the soft 
part of the dipoles: 



-r^soft 2)^13 -r^soft 2^13 

^12 3 = 1 \ N ' ^32 1 = -( \ N • (17) 

y\2{y\i-^yn)s\23 y23\y\i^yn)s\23 



Dokshitzer and Marchesini actually use slightly different dipoles, obtained upon averaging over 
the azimuthal angle of the emitted particle around the emitter. The differences in the choice for 
the dipoles will not be relevant to the rest of the paper. Note that although we label these terms 
"soft", the terms £>i2^3 and D^^ii have a soft and coUinear singularity. It is also clear, that in the 
coUinear limit does not reduce to Pq-^qg, but to 

(1-Z>12 ^ ^ 

In statements about shower algorithms reproducing the correct DGLAP equation in the coUinear 
limit, the term "correct DGLAP equation" refers therefore for the model above to eq. (fTSl) and 
not to eq. Q. 



3 Parton shower 

Whereas the factorization formula eq. ([7]) is exact in all singular limits, the r.h.s of eq. © does 
however not necessarily equal the full matrix element on the l.h.s away from the singular limit. 
(In the particular example discussed above the full matrix element for y* — > qgq actually equals 
the factorized form over the complete phase space, but this is not the general case.) Since for high 
parton multiplicities the full matrix elements are too complicated one approximates in a parton 
shower the full matrix elements by the factorized form over the complete phase space. Let us 
stress that this identification is exact in all singular limits and an approximation away from the 
singular limits. 

The antenna function A3 and the dipoles D123 and 1)32.1 are positive definite over the com- 
plete phase space and therefore can be interpreted as a probability distribution for the emission 
of an additional particle. For a parton shower algorithm we introduce two variables t ("shower 
time") and z ("momentum fraction"). The shower time t gives the scale at which the next split- 
ting occurs, the variable z describes for a splitting a ^ be the momentum fraction of the daughter 
b with respect to the mother a. For a complete description of a splitting we need in principle a 
third variable (|), but this variable will not be relevant for the discussion of this paper and we sup- 
press it. Therefore we can restrict ourselves to a two-dimensional space parametrized by (?,z) or 
(^12,^23) as in fig.[I] The choice for t (and z) is not unique. If we focus on the dipole £>i2,3 with 
singularities for ^'12 possible choices are 

^13 

for a I -ordered shower or 



■ln^ = -lnf^^^^Sl (19) 



^13 = -li^7i = -lM^i27g) (20) 



■^12 1 / ■^13 

for a virtuality-ordered shower. The quantity in eq. (fT9l ) and eq. (|20l ) is a fixed reference scale, 
usually taken to be the centre-of-mass energy squared of the showering system. The shower time 
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Figure 2: Lines of constant shower time t for a fc^-ordered shower (left) and a virtuality-ordered 
shower (right). The singular region for the dipole £)i2,3 is shown in blue and red. The singular 
region is approached for ? ^ oo. 

t takes values between a starting time tQ and +00. As a larger values of t corresponds to lower 
scales we have to require that the singular region is contained in the region defined hy t ^ 0°. In 
fig- m we show in the Dalitz plot lines of constant t for the definitions as in eq. ( fT9l) and eq. (|20|) . 
It is clear that the coUinear singular region yu = is contained for both definitions in the region 
defined by ? ^ 00. 

Associated to the scale t is the Sudakov factor Ai2,3, giving the probability that no emission 
occurs between the scales Iq and t 

Al2,3(?0,0 = j^dt' j d<Sfunres^{t' -ti^%nasCFDn3^ , (21) 

with given by eq. (fT9l) or eq. (I20l) . The derivative of the Sudakov factor with respect to t gives 
the probability of a splitting at the shower time t: 

d (X- f 

-^Ai2,3(?o,0 = l^^Cf d(^unres^{t-ti^) S123D 12,3 (22) 

Working this out for a virtuality-ordered shower and setting z = 1 — j23 one obtains 

'^Au.3M = ^Cf fdz(-^^-\-z]. d = ^e-'. (23) 



dt ' ' ' ' 271 7 Vl-z + 5 

5 



s 
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At finite shower time t the splitting probability is finite, as it should be. In the limit t ^ °o one 
recovers the DGLAP equation: 
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Figure 3: Lines of constant shower time t for an energy-ordered shower (left) and an angular- 
ordered shower (right). The singular region for the dipole £)i2,3 is shown in blue and red. For the 
energy-ordered shower lines of constant t intersect the singular region for finite t. For the angular- 
ordered shower the two singular points (^12,^23) = (O5O) and (^12,^23) = (O7 1) arc reached for 
finite shower time. 

A similar analysis can be carried out for a -ordered shower. 

In their paper Dokshitzer and Marchesini [35] did not use the transverse momentum or the 
virtuality as evolution variable. Instead they chose energy of the emitted gluon (in the rest frame 
of the dipole) as evolution variable. In our notation this amounts to the choice 

m - -ln|^-ln(^^f). (25) 

For this choice the lines of constant shower time t are shown in fig. [3l In this case t ^ 00 
corresponds to the single point yn = J23 = in phase space and lines of constant t intersect the 
coUinear region for finite t. We can therefore not expect to obtain a shower algorithm which is 
correct in the coUinear limit based on this evolution variable. Indeed, a short calculation shows 
in this case 




There are several problems related to an energy-ordered shower: For finite shower time t the 
splitting probability is infinite, due to the 1 / (ymax — yis) singularity in the integrand. In the limit 
t the integration over J23 reduces to a point and not to the integral over the splitting function. 
The integrand of eq. (|26l ) bears no resemblance to the Altarelli-Parisi splitting function. These 
deficiencies are all related to the inappropriate choice of the shower evolution variable. 
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For completeness we also show the corresponding plot for an angular-ordered shower in fig. [31 
Lines of constant angle are given by 

l-cosei2 = 2-. r (27) 

and the corresponding definition of the shower time t is given by 

= T^V (28) 

V(jl2+y23)(l-j23)eV 

The points (^12,^23) = (O5O) and (^12,^23) = (0, 1) are also reached for finite shower time t and 
angular-ordered showers have to introduce a cut-off on the variable z to avoid these points. In 
more detail we have for an angular-ordered shower 



-^Ai2,3(^0,0 = ^CfB{1 - ymax) J dy23 



1 







^ 1-^23 1 - ymax ^ J23 



3^23 1 - ( 1 - y23)ymax 1 " ( 1 " J23 )y« 



Q -t 



ymax = —e-'. (29) 
■^13 

As already mentioned above, we have for an angular ordered shower for finite t an infinite split- 
ting probability due to the soft singularity at 3^23 = 0- In angular-ordered shower programs this 
situation is usually handled by introducing an ad-hoc cut-off on the variable 3^23 ■ In the limit 
^ ^ 00 and setting z = 1 — J23 we recover the DGLAP equation: 

1 

liml--A,,3(,„,0l = S^-/*(t^->-^1- (30) 



4 The non-singlet quark fragmentation function 

In this section we review the relevant formulae for the evolution of the non-singlet quark frag- 
mentation function. The DGLAP evolution equation for the non-singlet quark fragmentation 
function d{x, Q^) is 

q'JqA^.q') = -/fic^^..W^(^'e2)' (31) 

X 

where Pqq is the regularized splitting function 

+ ^5(l-z)-(l+z). (32) 
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To a function f{x) we denote the Mellin transform by 

1 

f{N) = jdx^-'f{x). (33) 



In Mellin space the evolution equation for the non-singlet quark fragmentation function fac- 
torises: 



oc 

^'je2^>'2') = -:^CFPqq{N)d{N,Q^), (34) 

with 



Si (A'^ — 1) is the harmonic sum 



A'-l 2 



Si{N-l) = (36) 
,/=i J 

For large N, Si{N — I) diverges logarithmically. In Mellin space the evolution equation can be 
solved analytically 



d{N,Q') = (l + ^^Poln^) d{N^Ql), (37) 



where (3o = 1 1/3Ca — 4/3TrNf is the first coefficient of the beta-function. We are also interested 
in a toy model with = const, in this case the solution is given by 



d{N,Q') = (^^J d{N,Qt,). (38) 

The initial condition 

d{x,Ql) = 5(1 -x) (39) 

corresponds in Mellin space to 

d{N,Ql) = 1. (40) 

Finally we are interested in the .x:-space result for the quark fragmentation function d{x, (^) for 
values of x close to 1 and with = const. In this region only soft gluons have been emitted. 
With the ansatz [38] 

d{x.Q^) = A(e2)(i_^)«(e^) (41) 
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and neglecting terms which vanish in the limit x I one finds the equation 



1 dA{Q^) asCf 



+ ■ 



A{Q^) Jin 22 ;i 



--Y£-V|/(5(e2) + i) 



+ ln(l-jc) 



dB{Q^) a,CF 



0. (42) 



The coefficient of ln(l — jc) and the term independent of ln(l — x) have to vanish independently. 
With the initial condition eq. (|39l) one then obtains 



1 + 



In- 



21 



This solution is valid for 



ln(l-;c)- 



TT V 4 



Inrr-^lnCl. (43) 



'l-jc)<l and — In^— <1. 

71 1 — a: 



(44) 



5 Numerical studies 

5.1 The quark fragmentation function 

In this section we study the quark energy distribution in Monte Carlo events obtained from a 
shower simulation, starting from the hard matrix element e+e^ qq. As centre-of-mass en- 



ergy we take Q = mz, unless indicated otherwise. The (A'' 
distribution at the scale Qj is just 



d{N,Q 



l)-th moment of the quark energy 



(45) 



Our main interest is the comparison between the numerical shower program and the analytical 
result from the DGLAP equation. For this comparison it is sufficient to consider a toy model with 
= const. We set = 0.1. For the numerical result we first generate quark-antiquark events 
according to the hard matrix element, then start the shower at a scale Qo and run the shower 
to the lower scale 2ir. For a ^^-ordered shower starting from a process at the centre-of-mass 
energy Q the upper limit on Qq is given by Qo.max = 2/2. We then calculate the energy fraction 
of the quark (additional quarks obtained from g ^ qq splittings are not relevant to the discussion 
here): 



2E, 

Q 



(46) 



This defines a distribution in x. In addition, we can simultaneously bin the moments of this dis- 
tribution. 



We perform several comparisons between the numerical shower program and the analytical re- 
sult from the DGLAP equation. In Mellin space the DGLAP equation is an ordinary differential 
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Figure 4: Moments of the quark energy distribution at 2ir = 1 GeV obtained from starting tlie 
evolution at 2o = 2 GeV. 



equation and the numerical shower program has to reproduce this equation in the strongly or- 
dered limit. Strongly-ordered implies that the scale of successive emissions satisfy 

•••>e5-i>e5>e5+i>- m 

We can ensure these conditions by starting the shower at a rather low scale Qq = 2 GeV and run 
the shower only for a short interval to Qir = 1 GeV. The low starting scale Qo ensures ^ Qq 
and no emissions with a scale larger than Qq are generated. The short interval ensures that the 
number of events with two or more emissions is negligible and almost all events will have either 
zero or one shower emission. This is necessary since although in a shower successive emissions 
are ordered 

■■■>e--i>e->e-+i>-, (48) 

condition (1481) does not exclude successive emissions to be of the same order Qj = 0{Qj_i). 
Running the shower over the short interval from 2 GeV to 1 GeV tests therefore if the emission 
of a single particle correctly approaches the DGLAP limit. Fig. |4] shows for this case the com- 
parison between the numerical shower program and the analytical solution eq. (l38l) . We observe 
an excellent agreement. 
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Figure 5: Moments of the quark energy distribution at 2ir = 1 GeV for various centre-of-mass 
energies: Q = mz, Q = ITeV, Q = lOTeV and Q = lOOTeV. 



We then compare the numerical shower program and the analytical result from the DGLAP 
equation for a larger interval for the evolution. We start the shower at the hard scale Qq = Qo.max 
and run to the low scale Qir = 1 GeV. We compare again the moments of the quark energy 
distribution. We do this for the centre-of-mass energies Q = mz, 2=1 TeV, 2=10 TeV and 
2 = 100 TeV. We do not expect perfect agreement, since now the shower may generate emis- 
sions with a scale smaller but comparable to the previous one. However both the shower and the 
DGLAP equation resum the leading logarithm. In the limit where this logarithm is large against 
other terms, the results should agree. Fig.[5]shows the comparison between the numerical shower 
program and the analytical solution eq. (l38l) for the centre-of-mass energies Q = mz, 2=1 TeV, 
2 = 10 TeV and Q — 100 TeV. The starting scale of the shower is always 2o = Qo,max = 2/2. 
For the final scale of the shower the value 2ir = 1 GeV is always used. We observe that for large 
values of 

In^ (49) 

the two results approach each other. 
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Figure 6: The 10-th moment of the quark energy distribution for various values of a^^, obtained 
from a process with centre-of-mass energy Q = niz. The left figure shows the moment at Qir = 
1 GeV, the right figure shows the ratio between the values at Qir = 1 GeV and Qir = 2 GeV. 

Fig. [6] shows the 10-th moment of the quark energy distribution as a function of the (fixed) 
value of tty. The centre-of-mass energy is Q = mz and the starting scale of the shower Qq = 
Qo.max = 2/2. As the low scale Qir = 1 GeV is used. From eq. (l38l) we expect on a logarithmic 
scale a linear relationship with respect to the variation of a^: 



.|2C.^,,Wln|. 



(50) 



We observe in the left plot of fig. [6] that both the numerical result from the shower program and 
the theoretical curve give straight lines. However the slope is slightly different. From eq. (|50l) 
we see that the slope depends on the value of the hard scale Qq. To eliminate the dependence on 
Qo we show in the right plot of fig.[6]the ratio d{N, Qjji)/d{N,4Qjj^) of the 10-th moment at the 
low scales Qm and 22ir. In this ratio the dependence on Qq drops out: 



In 



fR) 



d{N,Q 
d{NAQ\^) 



—CFpqq{N)\nl. 



(51) 



We observe an excellent agreement. 



As a further comparison we now study the quark energy distribution in x-space for values of 
X close to 1. This region is sensitive to the emission of soft gluons. We start the shower at the 
hard scale Qq = Qo^max = 2/2 = '«z/2, and use as the final scale of the shower 2ir = 1 GeV. 
Fig. |7] shows the comparison between the numerical shower program and the analytical solution 
eq. (l43l) . We observe a good agreement. We would like to make a comment: The analytical 
solution gives a linear relation 



\nd{x,Q^) = A\n{l-x)+B, 



(52) 
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Figure 7: The quark energy distribution for small values of (1 
centre-of-mass energy Q = itiz- 



x) at 2iR = 1 GeV for the 



shown as a straight line in fig. Ul The validity of the analytical solution according to eq. (|44)) is 
restricted to the region 



- — <ln(l 
a. 



<0. 



(53) 



Similar, the shower does not generate emissions below a certain value of ln(l — Xmax)- Values 
of ln(l — x) below that value would correspond to emissions with a scale less than 2ir. As a 
consequence there are for any finite value of Qir events, which didn't radiate at all. (The fraction 
of these events is determined by the Sudakov factor at the scale Qir.) In fig.|7]we have normalized 
the shower result to the number of events which emitted at least one additional parton. 

To further test the validity of our conclusions under systematic variations of the shower as- 
sumptions, we use the capabilities of the ViNCiA plug-in to the Pythia 8 generator [5,22], 
which offers the possibility to use arbitrary radiation functions in conjunction with several differ- 
ent evolution variables within a dipole-antenna shower context. We also compare to the standard 
Pythia 8 p^-ordered shower [39], which represents a hybrid between the parton and dipole 
approaches. In both programs, we switch off g qq branchings, use a fixed ds = 0.1 and take 
the starting scale to be 2 = 10^ GeV as default, with a hadronization scale (infrared cutoff) of 
2iR = 1 GeV. In ViNCiA, we further switch off matching to matrix elements and matching to 
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Figure 8: Lines of constant shower time t for a -ordered shower (left) and an antenna-mass- 
ordered shower (right). The singular region for the antenna A (y 12, 3^23) is shown in blue and red. 
The singular region is approached for t ^ 00, 

2 — > 4 antenna functions. 

In ViNCiA, we consider two different dipole-antenna evolution variables: transverse momen- 
tum (identical to the Ariadne evolution variable [18]) and antenna-mass, referred to as Type I 
and Type II evolution, respectively. In terms of colour-ordered triplets of parton momenta, these 
variables are defined as follows: 

QKpuP2,P3)=^'-^^^pI. (54) 

5123 

2ii(pi,P2,P3) =2min(5i2,523) = '^ml^v (55) 

where the normalizations are chosen such that the maximal value of the evolution variable is 5'i23 
in both cases. Contours for constant values of these variables are shown in fig. [8l 

The most general form for a leading-log antenna function (dipole-antenna splitting function) 
for massless parton splitting is represented by a double Laurent series in the two branching 
invariants [24], 

AndsC ^ nt 3 ■^;7 

A(jl2,y23;*123) = — 2- 3^12 3^23 ' With 3;// = <1 , (56) 

*123 „^p^_i 5123 

where C is the colour factor. We here consider 3 different choices for these functions, the 
Gehrmann-de-Ridder-Glover ones ("GGG") [33], which are the defaults in ViNCiA, and system- 
atically high and low variations, MIN and MAX, respectively, with coefficients for the qq qgq 
and qg qgg functions given in Tab. [TJ We note that the default GGG qq qgq function is 
identical to the corresponding function in Ariadne and reproduces the Z qgq tree-level ma- 
trix element exactly. Also note that the colour factor for qg qgg is ambiguous up to l/N^ and 
that this variation is included in the MIN/MAX variation. 

In Pythia 8, the transverse-momentum variable, which we shall here call qr to distinguish 
it from the other definitions, agrees with the ARIADNE definition in the infrared limit, but differs 
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Figure 9: The quark fragmentation function as obtained with ViNCiA and Pythia 8, compared 
to the analytical expression, eq. (l43l) . Top Left: variation of the evolution variable (p^ -ordering 
vs. antenna- mass-ordering). Top Right: variation of the antenna functions (the Gehrmann- Glover 
functions vs. the ViNCiA MIN and MAX variations). Bottom Left: variation of the generator 
and of the hadronization cutoff contour (Pythia 8 vs. -ordered Vincia with a cutoff in p± 
or antenna-mass). Bottom Right: variation of the starting scale (mz = 10" GeV vs. niz — 10"^ 
GeV). In all cases, g ^ qq branchings were switched off, a constant = 0.1 was used, and D 
was normalized to the number of inclusive 3-parton events at the cutoff. 
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Table 1: Colour factors C (in the ViNCiA normalization [22]) and Laurent coefficients p for 
the antennae used in this study. The coefficients with at least one negative index are univer- 
sal (apart from a re-parametrization ambiguity for gluons). The positive-index coefficients are 
arbitrary and are here varied between MIN and MAX. 



from it by up to a factor of 2 away from that limit, see [39]. The splitting functions in Pythia 8 
are the ordinary DGLAP ones, augmented by matching to the tree-level Z qgq matrix element. 

In this study, we do not include variations of the kinematics maps beyond that offered by 
the default Pythia 8 and Vincia choices. Pythia 8 uses a partitioned-dipole-like map in 
which a "recoiler" recoils longitudinally (in the dipole centre-of-mass frame) against a "radiator". 
Vincia by default uses the gg ggg dipole- antenna map of Ariadne for all branchings. In 
the dipole-antenna case, no special distinction is made between radiator and recoiler; instead the 
proper coUinear limiting behaviour is obtained by a rotation angle going to 0° or 180° in the 
respective limits. 

In fig.[9l we show 4 plots illustrating the quark fragmentation function with these variations, 
as compared to the analytic expression, eq. (|43l) . The top left plot illustrates that one obtains good 
agreement between ViNCiA and the analytic expression, irrespective of the choice of evolution 
variable. (Recall that the validity of the analytic expression is limited to the region ln( l—x)<^ 
0.) It should be noted, though, that it could still be possible to obtain bigger variations for 
more extreme variations of the choice of evolution variable. The top right plot illustrates the 
dependence on the choice of antenna functions. As expected, this variation is larger, since this 
dependence enters already at first order in the shower expansion. Nonetheless, the asymptotic 
slope of all the curves agrees with the analytic expectation. In the lower left-hand plot, we show 
the variation between Pythia 8 and Vincia, with two different choices for the hadronization 
cutoff for the latter, either in or in antenna mass. Though we emphasize that this distribution 
is very infrared sensitive, the dependence on the choice of hadronization cutoff here seems rather 
mild. The PYTHIA 8 curve has a slightly different slope than the ViNCiA ones, but still appears 
to be within the uncertainty spanned by the variations above. Finally, in the lower right-hand 
plot, we show the results for two alternative starting scales, one at 100 GeV and the other at 
10000 GeV. As expected, the agreement improves with increasing energy (or, more precisely. 
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with increasing Q^j Qir). 

In summary, we find no evidence in either partitioned-dipole or dipole- antenna showers of 
a breakdown of agreement with DGLAP-based predictions of the quark fragmentation function, 
provided that "infrared sensible" evolution variables are chosen. The definition of infrared sen- 
sible is that both infinitely soft and coUinear emissions should be classified as unresolved for any 
finite value of the evolution variable. 



5.2 Comparison to second-order QCD 

A complementary check on the accuracy of the shower can be obtained by comparing its second- 
order expansion to second-order QCD matrix elements. This is simplest for the dipole- antenna 
shower, for which the number of possible histories for each phase space point grows less fast than 
for the partitioned-dipole case, so in this subsection, we shall only use dipole- antenna showers 
for the comparisons, but we emphasize that the results should be qualitatively similar for the 
partitioned-dipole case. This part of our study is similar to a previous comparison of Ariadne 
to second-order QCD by Andersson et al. [40]. 

In phase space regions dominated by leading logs, the ratio shower/matrix-element should be 
unity. In phase space regions dominated by hard wide-angle emissions, the shower could in prin- 
ciple be arbitrarily far from the matrix element, and finally in regions dominated by subleading 
logs (such as regions with two emissions at the same scale), the subleading-log properties of the 
shower can be probed. 

To perform this test independently of the shower generator, we use Rambo to generate a 
large number of evenly distributed 4-parton phase space points. For each phase space point, we 
evaluate the leading-colour 4-parton antenna function, 

^^^^ ^ — wm — ' ^ ^ 

as given by Gehrmann et al. [33] (counter checked with [41] to protect against typos). 

We then compute the tree-level leading-colour LL antenna- shower approximation corre- 
sponding to the same phase space point, based on nested 2^3 branchings. For 4 partons, 
there are two possible antenna- shower histories; 

• A) parton 2 emitted between partons 1 and 3. 

The 4-parton evolution scale is then = 2|(1,2,3). 

• B) parton 3 emitted between partons 2 and 4. 

The 4-parton evolution scale is then 2|g = 2|(2, 3, 4), 

with Qe denoting a generic evolution variable. We shall here consider energy-ordering, px.- 
ordering, and antenna-mass ordering. 

We note that a similar study for parton- or partitioned-dipole (Catani-Seymour) showers 
would need to consider 8 possible paths from 2 to 4 partons [42]: two possible radiators in 
the first 2 — > 3 step (the quark and the antiquark), and 4 possible radiators in the subsequent 
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3^4 step (treating the two "sides" of the gluon, which are generally associated with differ- 
ent kinematics mappings, as separate). The dipole- antenna shower is thus very economic in the 
number of terms generated at each successive order. 

Using a new clustering algorithm that contains the inverse of the ViNCiA 2^3 kinematics 
maps, we may perform clusterings of the type (a, r,b) — > {d,b) in a way that exactly reconstructs 
the intermediate 3-parton configurations that would have been part of the shower history for each 
4-parton test configuration, for each of the paths A and B0. This gives us an exact tree-level 
reconstruction of how the antenna shower would have populated each path. We can now use this 
to test the shower approximation over all of 4-parton phase space. 

We shall do this by plotting the ratio 



(12,23,4)A,g(l,2,3)+A,^(l,23,34)A,^(2,3,4) 
A4(l, 2,3,4) 



with hatted variables ij denoting clustered momenta. R thus gives a direct measure of the amount 
of over- or under-counting by the shower approximation, with values greater than unity corre- 
sponding to over-counting and vice versa. 

The ratio above, eq. (|58]) . contains nested products of antennae, identical to the subtraction 
terms that would be used in a fixed-order calculation. This does not take into account the ordering 
condition imposed in the shower, however. To impose this condition, we must include step 
functions in the shower approximation, as follows 

„£ 0(23J-e«)A,^(12,S,4)A,j,(l,2,3)+0(e3jj-e4B)A,^(l,S,34Mj,^(2,3,4) 

= A4(l,2,3.4) • 

where 

Q4A = e£(l,2,3) ; Q3A = e£(12,23,4) 
Q4B = 2^(2,3,4) ; Q3B = ^£(1, 23,34) 

The ratio Rf now faithfully reproduces the shower approximation expanded to tree level, phase 
space point by phase space point, for an arbitrary choice of evolution variable, Qe- 

Since the full 4-parton phase space has more dimensions than can fit on paper, and since the 
leading singularity of the gluon emission antenna functions goes like (with defined as 
in eq. (|54l) ). we project the full phase space onto two p± values, one of which we choose to 
correspond to the initial 2 — 3 step of a would-be shower history and the second to the 3 ^ 4 
step. Specifically, the ordinate along the y axis will be 

3; axis: j3i;2 = min(p^(l,2,3),pi(2,3,4)) , (61) 

corresponding to the second branching, and the ordinate along the x axis will be the p^ value of 
the reclustered 3-parton configuration corresponding to the min(...), 

._;i2,23,4) ;pi(l,2,3)< Pi (2,3,4) 
xaxis: P±i = \ ■ (62) 

;i, 23,34) ;pi(2,3,4)<pi(l,2,3) 




'Note that the inversion of ViNCiA by this clustering algorithm is exactly one-to-one, with no approximation 
made. This was validated by reclustering a large number of actual branchings generated by the shower and recover- 
ing the pre-branching configurations exactly, including global orientations, etc. 
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Each point in (/?i;i,/'±;2)-space thus contains a distribution of all 4-parton configurations 
with that particular combination of p±-j and p±-^2 values. We shall plot both the average of this 
distribution, which we call < >, as well as a measure of the spread of the distribution, which 
we define as 

RMS(7?f ) = ioRMS(iogio(4)) > 1 . (63) 

By using this form we probe the average factor of deviation from unity rather than the absolute 
measure of the deviation itself. (I.e., we want a point with Rf =0.1 to count as having a deviation 
of a factor of ten, rather than an absolute deviation of 0.9, from unity.) A special case is when 
we encounter dead zones in which the shower answer is zero, and hence the factor of deviation 
would nominally be infinite. When computing the RMS above we therefore put a floor on the 
deviation at a factor 0.01 times the matrix element. 

In figs. [To] and [TTJ we show the average and RMS for four different ordering variables; for 
comparison, we first show the result without any ordering imposed, as in eq. (I58l) . i.e., a simple 
product of nested antennae with no functions imposed. This is equivalent to the subtraction 
terms constructed for fixed-order calculations and can be represented by "ordering" in the vari- 
able m^ax = max (5 12, 523) • We then compare to the ordered results, eq. (|59| ). for energy-ordering 
(as defined by Dokshitzer and Marchesini, i.e., ordering in the energy of the emitted parton in the 
CM of the Z boson), p^-ordering (as defined in eq. (I54l)). and antenna-mass-ordering (as defined 
in eq. (|55])). We may identify several regions of interest on the plots shown in figs. [TOl and [TTl 

• Origo: double-LL singular region: -C p±-i <ti s, i.e., two widely separated jets plus 
two strongly ordered emissions. Should be correctly described by any LL shower. 

• Top right-hand corner: Hard region: p±-i ~ p±-2 ~ i-C, 4 widely separated jets. Should 
be correctly described only by the 4-parton matrix element. 

• Area below diagonal dashed (cyan) line (bottom right-hand comer): region in which the 
second emission is strongly ordered with respect to the first: p±-2 P±;i- Should be 
correctly described by any LL shower matched to the 3-parton matrix element. 

• Intersection of y-axis with diagonal solid (cyan) line: Single-NLL (double-emission) re- 
gion: p^.2 ~ p±-i <S s, i.e. 2 widely separated jets plus one strongly ordered 2^4 emis- 
sion (two powers of but only one large scale difference). Should be correctly described 
only by the 4-parton matrix element and/or by an NLL shower. 

• Area above diagonal solid (cyan) line (except top right-hand corner, see above): pj^-unor- 
dered region. Corresponds to a dead zone in a shower ordered in p±. Although this zone 
occupies a relatively large area in our projection, this is chiefly an artifact of our choice of 
variables. The actual phase space volume in this region amounts to roughly 1.5% of the 
full 4-parton phase space. 

Our operational definition of a "correctly described region" we shall here take to be that both 
the average of the Rf distribution as well as its RMS factor in that region should tend to unity. 
If so, this means that the shower is not only getting the average of the distribution right, but 
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Figure 10: Left: Average of the Rf distributions, eqs. (l58l) and (|59| ). for no ordering (top) and 
energy-ordering (bottom). Diagonal lines indicate boundaries between unordered, ordered, and 
strongly ordered regions (doubly strongly ordered region is at origo). Right: the RMS of the 
factor of deviation from unity of the Rf distributions. 
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Figure 11: Left: Average of the Rf distributions, eq. (|59| ). for p^-ordering (top) and antenna- 
mass -ordering (bottom). Diagonal lines indicate boundaries between unordered, ordered, and 
strongly ordered regions (doubly strongly ordered region is at origo). Right: the RMS of the 
factor of deviation from unity of the Rf distributions. The RMS distribution for antenna-mass- 
ordering (on the lower right) is somewhat affected by the occurrence of a few dead points in a 
region extending towards the lower left, which "artificially" increase the RMS in that region. 
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there are also no large fluctuations on either side of the average. If only the average is unity but 
the RMS factor is not, then the interpretation is that the shower is still over- and under-counting 
individual phase space points, and hence the relevant part of phase space is not being populated 
accurately. 

In the top row of Fig.[lOl we show the average and RMS factor without imposing any ordering 
condition at all (apart from that implied by the nested phase spaces). The central grey colour 
towards the bottom right indicates that the shower approximation deviates less than 10% from 
the matrix element there, and the two neighbouring grey shades indicate 20% deviation. The 
next colours (lighter red and darker blue shades) represent factors of 2, 5, and 10, respectively. 
The dominance of light shades in the upper left half of the plots thus indicate that the nested 
LL antenna products, without ordering, exhibit a large over-counting whenever the 2"'^ emission 
does not have a several orders of magnitude smaller than the one. 

As expected, the situation for energy-ordering is actually worse (bottom plot row). Though 
the agreement improves slightly for hard 2"'^ emissions, we see a disturbing trend towards sys- 
tematic undercounting in the double-LL region, indicated by the darkening shades towards origo. 
What is even more disturbing is that the RMS factor does not follow the average, but instead 
blows up towards the LL singular regions. This means that the "good" average agreement on 
the left-hand side is really obtained by cancellations between large over- and undercounting of 
individual phase space points. The desired LL singular behaviour is therefore observed not to be 
obtained for this evolution variable. 

Finally, the quite impressive properties of j!7^-ordered dipole- antenna showers are evident in 
the near-unity average value of both and its RMS factor of deviation over all LL-dominated 
phase space regions in the top left-hand plot of fig. El as was also noted in the study of Anders son 
et al. [40]. Though we here use a different set of antenna functions and kinematics maps, we see 
that the excellent basic agreement with the second-order QCD matrix elements is retrieved also 
in our case. We note that by definition a -ordered shower does not generate any points in the 
/(^-unordered region above the diagonal (cyan) solid line. (Again, these dead zones were also 
pointed out in the study by Andersson et al. [40]). As mentioned, we find that roughly 1.5% of 
the full 4-parton phase space is left unpopulated by this particular ordering variable (antenna- 
mass-ordering gives a similar number). Keep in mind that this is 4-parton space though; none of 
these showers have any dead zones in 3-parton phase space. 

The case of antenna-mass-ordering (fig. [TTl bottom row) is similar to jc^ -ordering, but its 
dead zones do not follow strict contours of p^, and hence the RMS factor looks "artificially" 
large over that part of phase space in which dead points exist, roughly the area above the diagonal 
of the plot. 

We plan to return to the issue of dead zones in a future paper, but note that since they are 
located in the NLL-dominated region, they do not affect the conclusions we wish to make here 
concerning the LL behaviour of the evolution choices. We thus restrict ourselves to the conclu- 
sion that pj^-ordered dipole- antenna showers appear to give an excellent approximation to the 
full 4-parton matrix element over all LL dominated regions of phase space (below the diagonal 
(cyan) dashed line). When -ordering is imposed, the RMS factor of deviation furthermore 
registers an impressive sharpening-up of the i?f distribution, yielding much larger regions of 
unity RMS factor than the corresponding case without ordering. Also antenna-mass-ordering 
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represents a substantial improvement, although the improvement in the RMS is slightly masked 
by the fact that our projection "smears out" its dead zones over a larger area of the plot than for 
-ordering. 

Energy ordering, on the other hand, effectively introduces artificially under-counted zones 
in the doubly-LL singular region, while still not removing the over-counting that was already 
present in the same region without ordering — hence the RMS measure of deviation actually 
worsens as we go further into the singular region. The factorization implied by this choice of 
evolution variable is thus clearly not consistent with the structure of QCD. 

6 Conclusions 

In this paper we studied shower algorithms based on partitioned-dipoles and dipole-antennae. In 
particular we investigated the behaviour in the coUinear limit and showed that with an "infrared- 
sensible" definition of the evolution variable they reproduce the DGLAP evolution equation. The 
definition of infrared sensible is that both infinitely soft and coUinear emissions should be clas- 
sified as unresolved for any finite value of the evolution variable. Examples of such choices are 
fcj^-ordering or mass-ordering (ordering in virtuality for partitioned-dipole and antenna-mass for 
dipole-antenna showers, respectively). On the other hand, ordering in the energy of the emitted 
particle is not infrared sensible (it classifies infinitely coUinear emissions as being resolved) and 
does not reproduce the DGLAP equation. 

In addition to these analytic arguments, we have also presented a numerical study, making use 
of existing dipole shower algorithms. We demonstrated that the DGLAP behaviour of the quark 
fragmentation function is reproduced by these models for a range of different infrared sensible 
shower algorithms, in particular p^-ordered ones. In addition we compared dipole-antennae to 
second-order QCD matrix elements and again retrieve good agreement in the strongly ordered 
(LL-dominated) region for -ordering, but not for energy-ordering. 
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